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O ! Abstract 

(N ■ 

Oh 1 Extremely asymmetrical scattering (EAS) is an unusual type of Bragg scattering 

rj^ ■ in slanted periodic gratings with the scattered wave (the +1 diffracted order) 

propagating parallel to the grating boundaries. Here, a unique and strong sen- 
sitivity of EAS to small stepwise variations of mean structural parameters at the 
1 • grating boundaries is predicted theoretically (by means of approximate and rig- 

orous analyses) for bulk TE electromagnetic waves and slab optical modes of 
arbitrary polarization in holographic (for bulk waves) and corrugation (for slab 
q ■ modes) gratings. The predicted effects are explained using one of the main phys- 

ical reasons for EAS — the diffractional divergence of the scattered wave (similar 
to divergence of a laser beam). The approximate method of analysis is based on 
this understanding of the role of the divergence of the scattered wave, while the 
rigorous analysis uses the enhanced T-matrix algorithm. The effect of small and 
Qh| large stepwise variations of the mean permittivity at the grating boundaries is 

analysed. Two distinctly different and unusual patterns of EAS are predicted in 
the cases of wide and narrow (compared to a critical width) gratings. Compari- 
^ | son between the approximate and rigorous theories is carried out. 

-xt- : 

0\ '. 1 Introduction 

O ; 

Previously, it has been demonstrated that scattering in slanted, strip-like, wide (compared to the 
wavelength) periodic gratings with the scattered wave propagating parallel or almost parallel to the 
front grating boundary is characterized by a strong resonant increase in the scattered wave amplitude 
[1-12]. This type of scattering was called extremely asymmetrical scattering (EAS). It has been shown 
to be radically different from the conventional Bragg scattering in transmitting or reflecting gratings 
[1-12]. For example, in addition to the strong resonant increase of the scattered wave amplitude at the 
Bragg frequency, scattering in the geometry of EAS may also be characterized by additional unique 
strong resonances, such as an additional exceptionally strong resonance in the sidelobe structure at 
a frequency that is noticeably higher than the Bragg frequency [11], additional unique resonance 
with respect to angle of scattering if the scattered wave propagates almost parallel to the grating 
boundaries [12], combinations of strong simultaneous resonances in non-uniform gratings [7-9], etc. 
It has also been shown that the diffractional divergence of the scattered wave inside and outside the 
grating (similar to the divergence of a laser beam of finite aperture) is responsible for all these unique 
wave effects in the geometry of EAS [2-9,11,12]. The necessity of taking the divergence into account 
can be seen from the fact that if the divergence is neglected and the scattered wave propagates parallel 
to a strip-like grating (the geometry of EAS), then this wave must be located within this grating. Thus 
we would have a scattered beam with aperture equal to the grating width even if the incident wave is 
a plane wave. This scattered beam must spread outside the grating owing to diffractional divergence. 
For a more detailed analysis of the role of diffractional divergence in EAS see [2-9,11,12]. Note that 
when the scattered wave propagates at a significant (usually several degrees) angle with respect to 



the grating boundaries, the diffractional divergence of the scattered wave can be neglected [12], and 
we have conventional scattering in transmitting or reflecting gratings. 

On the basis of understanding the role of the diffractional divergence for scattering in the EAS ge- 
ometry, a new powerful method of simple analytical (approximate) analysis of this type of scattering 
has been introduced and justified [2-9,11,12]. The main advantage of this method is that it is di- 
rectly applicable to the analysis of all types of waves in different periodic gratings, including surface 
and guided waves in periodic groove arrays [2-9,11,12]. Even in a complex five-layer structure with 
non-collinear grating-assisted coupling, the new method allowed accurate analytical analysis of ex- 
tremely asymmetrical coupling of two optical modes guided by neighbouring optical slabs with a 
corrugated interface [2]. The analysis of the applicability conditions for the new method, based on 
physical speculations [12] and rigorous analysis of EAS [10], has demonstrated that this method nor- 
mally gives excellent agreement with the rigorous theory for gratings with small amplitude. It is 
these gratings with small amplitude that are of most interest from the viewpoint of EAS since they 
result in strong resonant increase of the scattered wave amplitude [2-9]. In addition, the new method 
provides excellent insight into the physical reasons for EAS, which will allow thoughtful selection of 
optimal structural parameters for future EAS-based devices and techniques. 

If required (for example, for very narrow gratings or large grating amplitude), rigorous analysis of 
EAS of bulk electromagnetic waves can also be carried out by means of one of the known numerically 
stable rigorous approaches, such as an enhanced T-matrix approach [13,14], S-matrix and R-matrix 
approaches [15,16], an S-matrix approach using an oblique Cartesian system of coordinates [17], a 
C method [18] (for EAS of guided modes in corrugation gratings), etc. For example, in a recent 
paper [10], the enhanced T-matrix approach was used for the rigorous analysis of EAS of bulk TE 
electromagnetic waves in uniform holographic gratings with constant mean dielectric permittivity. 
Note, however, that all these rigorous numerical methods do not reveal the main physical reason for 
the unique pattern of scattering in EAS — the diffractional divergence of the scattered wave. 

Since the divergence plays a crucial role for EAS, variations in this divergence are expected to result 
in substantial variations in the whole pattern of scattering. It is well known that stepwise or gradual 
variations of mean structural parameters across a laser beam (for example, in nonlinear wave prop- 
agation [19]) may result in substantial variations in the diffractional divergence. Therefore, it can be 
expected that EAS will be unusually sensitive to even small variations of mean structural parameters 
across the grating. 

On the one hand, this unusual sensitivity may present a problem for experimental observations of 
EAS, since variations of mean structural parameters naturally occur during manufacturing of peri- 
odic gratings. For example, etching or ruling processes that are used for the fabrication of strip-like 
relief gratings on a surface of a planar waveguide usually result in a small reduction of the mean 
thickness of the waveguide (slab) in the region of the grating. Thus, mean thickness of the guid- 
ing slab will experience a step-like variation at the front and rear grating boundaries. Fabrication 
of holographic gratings for bulk and guided optical waves also results in varying mean dielectric 
permittivity in the grating. For example, if coherent UV radiation is used for writing a grating in a 
photosensitive material, the mean dielectric permittivity in the region of the grating will be slightly 
increased compared with the regions outside the grating that are not affected by the UV radiation. 
Such small variations of mean structural parameters can normally be ignored in the case of conven- 
tional Bragg scattering in transmission and reflection gratings. However, in the case of EAS, these 
variations may result in very significant changes in the pattern of scattering. 

On the other hand, the unusual sensitivity of EAS to small variations of mean structural parame- 
ters may be very useful, e.g. for the development of new optical and acoustic sensors and precise 
measurement techniques. For example, consider EAS of an optical slab mode into another guided 
mode of the same slab in a strip-like grating on one of the slab interfaces. Deposition of thin films 
or layers of some substance onto the regions of the slab outside or inside the grating will result in 
varying mean effective permittivity of the slab at the grating boundaries. Thus, EAS can be used for 
the detection of such deposited layers (thin film sensors). 

However, no theoretical analysis of EAS of optical waves in non-uniform gratings with varying mean 



structural parameters has been carried out so far. 

Therefore, the aim of this paper is to present theoretical (approximate) and numerical (rigorous) anal- 
yses of EAS of optical bulk and guided waves in gratings with step-like variations of mean dielectric 
permittivity (for bulk TE waves) and mean slab thickness (for optical guided TE and TM modes) at 
the front and rear grating boundaries. The theoretical analysis will be based on the recently devel- 
oped approximate method allowing for the diffractional divergence of the scattered wave [2-9,12]. It 
will be used for the investigation of EAS of bulk and guided optical waves in holographic and surface 
relief gratings, respectively. The rigorous numerical analysis will be presented for bulk TE electro- 
magnetic waves and will be based on the enhanced T-matrix approach developed by Moharam et al. 
[13,14]. Comparison between the approximate and rigorous results will be carried out. In particular, 
the unusually strong sensitivity of EAS to small variations in the mean structural parameters at the 
front and rear grating boundaries will be confirmed and analysed. 



Figure 1: The geometry of EAS in a slanted holographic grating of width L with stepwise variations 
of the mean permittivity at the front and rear boundaries. The vectors kn, k2i, and are the wave 
vectors of the scattered wave (the first diffracted order) in front, inside, and behind the grating. For 
EAS of slab modes in a corrugation grating, the plane of the figure is the plane of the slab, and £\, 
£2, £3, and e g must be replaced by the slab thicknesses h\, hi, the mean slab thickness in the grating 
region /z 3 , and the corrugation amplitude £, g/ respectively. 

The structure analysed in this paper is presented in figure 1. First consider bulk electromagnetic 
waves in a holographic grating with sinusoidal variations of the dielectric permittivity within a slab 
of thickness L: 



where the coordinate system is shown in figure 1, £2 is the mean dielectric permittivity in the grating, 
£\ and £3 are the dielectric permittivities of media 1 and 3 surrounding the grating (figure 1), £ g is the 
complex amplitude and q = (q x , q y ) is the reciprocal lattice vector of the grating that is parallel to the 
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Xo-axis, q = In I A, A is the period and L is the grating width. The grating is assumed to be infinite 
along the y- and z-axes. It is also assumed that dissipation is absent, i.e. e 1/2 ,3 are real and positive 
(the effect of dissipation on EAS is considered in [20]). The variations of the mean permittivity at the 
front and rear boundaries are given byAei = £\ — e 2 and Aeo, = £3 — e 2 (these values can obviously 
be real positive or negative). A TE electromagnetic wave with the amplitude of the electric field S10 
and wave vector ki is incident onto the grating at an angle Q\q in the xy plane — figure 1 (non-conical 
scattering). 

The solution to the wave equation in the grating can be written in the form of the expansion [13,14,21]: 

+00 

E 2 (x,y,t)= £ E 2n (x,y,t), (2) 

n— — 00 

where 

E 2 „(x, y, t) = S 2n exp(ixk 2nx + iyk 2ny - iwt) 

is the field in each of the inhomogeneous waves in the sum (2) with the x-dependent amplitude 
S 2n (x), and k 2nx and k 2ny are the components of the wave vectors 

k 2 „ = k 20 -nq, (n = 0, ±1, ±2, ...). (3) 

1 /2 

If for some value of n the magnitude of the wave vector k 2n is equal to tve 2 ' /c, where w is the 
frequency of the wave, and c is the light speed in vacuum, then the Bragg condition is satisfied for 
this value of n, and the nth wave in expansion (2) may have amplitude that is comparable with or 
even larger than the amplitude of the incident wave (zeroth term in equation (2)). 

In this paper it is assumed that the Bragg condition is satisfied precisely for the first diffracted order, 
i.e. for n = 1: 

k 2 i = k 2 o - q, 

1 /2 

where k 2 \ = |k 2 i \ = we 2 / c. The wave vector of the scattered wave (the first diffracted order) in the 
grating, k 2 i, is parallel to the grating boundaries, i.e. parallel to the y-axis (figure 1) — the geometry 
of EAS [1-10]. 

In the approximate method of analysis of EAS [2-9,11,12], only the zeroth order (incident wave) and 
the first order (scattered wave) in equation (2) are considered to be significantly non-zero (the two- 
wave approximation). It has been shown that the allowance for the second-order x-derivative of the 
scattered wave amplitude in the grating is essential for the correct description of scattering in the 
geometry of EAS [2-9,11,12]. In the approximate method [2-9,11,12], this second-order derivative is 
introduced into the coupled wave equations through consideration of the diffractional divergence of 
the scattered wave, which is one of the main physical reasons for resonant wave effects in the EAS 
geometry [2-12]. This approach is especially useful for the analysis of EAS of surface and guided 
waves [2-9, 12], where other methods fail to provide any reasonably simple solutions. 

The approximate method developed in [2-9,12] is also directly applicable to the case of non-uniform 
gratings with step-like variations of the mean structural parameters at the grating boundaries. For ex- 
ample the coupled wave equations in the grating in the geometry of EAS can be written as [4,6,9,12]: 

d 2 S 21 (x)/dx 2 + K S 20 (x) = 0, 
&S 20 (x)/dx-iK x S 2l (x) = 0, (4) 

where 

K = -2k 21 V sin(n - 2O ), K x = r\ cos(rj)/ cos)0 2O ), (5) 

9 2 q is the angle of refraction of the incident wave at the front boundary, 77 is the angle measured from 
the x -axis to the wave vector of the incident wave k 20 in the grating counter-clockwise (figure 1), T 
and Ti are the coupling coefficients that are determined in the conventional coupled wave theories 
for non-slanted gratings with fringes parallel to the grating boundaries. 

The incident wave is partly reflected from the front grating boundary, because the mean dielectric 
permittivity experiences a step-like variation at x = 0. However, in the considered approximation, 



this process of reflection (and refraction) is independent of scattering inside the grating. Therefore, 
the reflected and refracted waves at the front boundary x = can be found separately from scattering 
by means of the well-known Fresnel's equations for an interface. The determined amplitude S200 of 
the refracted wave at the grating boundary x = 0, and the angle of refraction 020 must then be used 
in equations (4) and (5) as the amplitude and the angle of incidence of the incident wave at this 
boundary. The wave vector of the incident wave in the grating should be taken as k 2 o — the wave 
vector of the refracted wave. 

Similarly the reflection of the incident wave at the rear boundary can also be considered separately 
from the process of scattering. In this case, the incident wave in the grating at x = L is regarded 
as the wave transmitted through the grating. Then the reflection from the rear boundary will be the 
reflection of the transmitted wave from this boundary which (in the approximate method of analysis 
[2-9,12]) has no effect on the process of scattering. 

It can also be seen that if £ 2 7^ £\, £3, then the scattered wave outside the grating should not neces- 
sarily propagate parallel to the grating boundaries, even if this is the case inside the grating (figure 
1). In general, the solutions for the scattered wave outside the grating are 

E\\(x,y r i) = Aq exp{\k\\ y y — ikllxx — icvt), 

E 31 (x,y,t) = B exp(ik 31y y + ik 31x x-iwt), (6) 
where A and Bo are the amplitudes of the scattered waves in media 1 and 3, respectively, k lly = 

k nx = (k 2 n - klrfl 1 , k 3lx = (k 2 31 - (7) 

1 /2 1 /2 

where k\\ = a>£ 1 / c, k 3 \ = a>£ 3 / c, and k\\ x and k 3 \ x are chosen positive real or positive imaginary. 
Note that in the rigorous analysis, equations (6) are replaced by the Rayleigh expansions of the elec- 
tric field in the regions x < and x > L [13-16,18,21]. In this case the waves with amplitudes Aq and 
Bo (equations (6)) are the first diffracted orders in these expansions. 

In the two-wave approximation, the boundary conditions at the grating boundaries can then be writ- 
ten as [2-9,12]: 

S2o\x=o = S2oo/S2iU=o = ^0/ (dE 2 i/dx) x= o = (dEu/dx) x=0 , 

E 2 i\ x =l = E 3 i\ x=L ,(dE 2 i/dx) x=L = (dE 31 /dx) x=L . (8) 

Note again that the refraction of the incident (transmitted) waves at the front and rear boundaries of 
the grating has not been considered when writing boundary conditions (8). 

For the rigorous theory, the boundary conditions at the grating boundaries are given in [13,14,21]. 

Boundary conditions (8) determine unknown wave amplitudes A 0/ B , and three constants of inte- 
gration in the solutions to the coupled wave equations (4) in the grating region [2-6]. These constants 
determine the incident and scattered wave amplitudes inside and outside the grating (for the explicit 
form of the solutions see [2-6]). 

3 Numerical results 

The numerical results of the approximate and rigorous analyses of EAS in gratings described by 
equation (1) are presented in this section for bulk TE electromagnetic waves. As has been indicated 
in [6-9,11,12], there are two typical patterns of EAS, which correspond to narrow and wide gratings. 
Narrow gratings are those whose widths are less than a critical width L C/ wide gratings are those 
with widths larger than L c [6-9,11,12]. Physically, L c /2 is the distance within which the scattered 
wave (beam) can be spread across the grating (i.e. in the direction normal to the wave propagation) 
by means of the diffractional divergence, before being re-scattered by the grating [7-9]. Simple meth- 
ods of determination of the critical width were developed in [7, 9]. This critical width is extremely 
important for scattering in the geometry of EAS. It appears to determine numerous resonant effects 
in uniform and non-uniform gratings [7- 9,11,12]. Strong differences in the patterns of scattering will 
also be demonstrated here for EAS in narrow and wide gratings with step-like variations of the mean 
structural parameters. 



3.1 Narrow gratings 
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Figure 2: The x-dependencies of relative amplitudes | S21 (x) / S10 1 of the scattered bulk TE electro- 
magnetic wave in the holographic grating with the parameters: £2 = 5, e g = 5 x 10~ 3 , L = 10 n, 
010 = 7r/4, A(vaccum) = l\i, the grating period (A ps 0.584^im) and the orientation of the 
fringes are determined by the Bragg condition and the direction of the wave vector k 2 i. (a, b) 
\Deltae\ = e\ — e 2 7^ and \Deltae3 = e 3 — e 2 = (variations of the mean permittivity only at 
the front boundary); (c, d) \Deltae\ = e\ — e 2 = and \ Deltae3 = e 3 — e 2 7^ (variations of the mean 
permittivity only at the rear boundary). Curves 1: Aej = and Ae 3 = (no variations of the mean 
permittivity). Curves 2, 3, 4, 5: (a) Ae t = -10 -5 , -10 -4 , -10 -3 , -KT 1 , respectively; (b) Ae x = 1(T 5 , 
lCT 4 , -CT 3 , 10" 1 ; (c) Ae 3 = -1(T 5 , -KT 4 , -1(T 3 , -KT 1 ; (d) Ae 3 = KT 5 , 1(T 4 , -CT 3 , 10" 1 . Curve 
6: Aei = —4, i.e. £\ = 1 (vacuum in front of the grating). 

Figure 2 presents the dependencies of the normalized scattered wave amplitude on the x-coordinate 
inside the holographic grating with L = lOjtmi, e 2 = 5, e g = 5 x 10~ 3 , A(vaccum) = ljum, d\o = 45°, 
and step-like variations of the mean permittivity at the front (figures 2(a), (b)) and rear (figures 2 
(c), (d)) grating boundaries as indicated in the figure caption. The Bragg condition is assumed to be 
satisfied precisely. The above structural parameters correspond to the critical width L c SOjUm [7,9], 
and therefore, the dependencies in figure 2 are typical for narrow gratings with L < L c . 
The most important feature that can be seen from this figure is that the scattered wave amplitude 
is unusually sensitive to small variations of the mean permittivity at the grating boundaries. Even 



very small step-lik e variations of the mean permittivity at either of the boundaries (Aei 3 = 10~ 5 = 
2 x 10~ 6 e 2 <C e g = 10~ 3 e 2 ) may result in noticeable changes in the scattered wave amplitude in 
the grating — compare curves 1 and 2 in figures 2(b), (d). Note that further decrease of grating width 
results in an approximately proportional increase of sensitivity of EAS to small variations of the mean 
dielectric permittivity at the grating boundaries. This is related to sharper EAS resonance in narrower 
gratings (if the width is less than L c ) [4,5,12]. 

The unusually strong sensitivity of EAS to small variations of mean structural parameters has no 
analogies in conventional Bragg scattering. As has been mentioned in section 1, it can be explained 
by a strong sensitivity of the diffractional divergence of the scattered wave to small variations of 
structural parameters along the wave front. 

Another interesting aspect is that varying mean permittivity at one of the grating boundaries results 
in significant changes of the scattered wave amplitude throughout the structure (figures 2(a)-(d)). 
Moreover, if the dielectric permittivity outside the grating is larger than inside (i.e. Aei,3 > — see 
figures 2(b), (d)), then the scattered wave amplitude is almost constant across the grating (especially 
for small values of Ae^). If Aei,3 < (figures 2(a), (c)), then the dependence of the scattered wave 
amplitudes on the x-coordinate is fairly noticeable. However, if the grating width is reduced further 
(e.g. down to « 5^xm), then the scattered wave amplitude becomes approximately constant across 
the grating for both the cases with Aei 3 > and Aei 3 < 0. This is because in narrow gratings 
(with L < L c ) the diffractional divergence is very efficient in spreading the scattered wave across the 
grating (see also [6-9]). Thus, any changes in the scattered wave (caused, for example, by varying 
mean structural parameters) at one boundary of the grating are effectively felt at the other boundary 
(and everywhere in the grating). Therefore, variations of the scattered wave amplitude in the grating 
are smoothed out by the diffractional divergence, and this results in only weak dependence of the 
scattered wave amplitude on the x-coordinate (see curves 1-4 in figures 2(a)-(d)). If the grating width 
is decreased, the diffractional divergence becomes even more efficient (at shorter distances), and the 
x-dependencies of the scattered wave amplitude become even weaker for all values of Aei 3. 
It can also be seen that if Aei ,3 > (i-e. the dielectric permittivity outside the grating is larger than 
inside), then the sensitivity of EAS to small variations of the mean permittivity at the grating bound- 
aries is stronger than for Aei ,3 < (compare curves 1-4 in figures 2(a), (c) with curves 1-4 in figures 
2(b), (d)). This is because if Ae^ < 0, then the waves outside the grating are exponentially decaying 
with increasing distance from the grating (see also section 2). These waves cannot be associated with 
an energy flow towards or away from the grating. On the other hand, if Ae^ > (figures 2(b), (d)), 
then the scattered waves outside the grating (according to Snell's law) are propagating waves travel- 
ling away from the grating. The resultant energy losses from the grating cause more significant (than 
in the case of Aei^ < 0) reduction in the scattered wave amplitudes (figures 2(a)-(d)). The smaller 
the grating width and /or grating amplitude e g , the stronger the difference in sensitivity of EAS to 
small positive and small negative values of Aei^. 

Note that the approximate and rigorous analyses of EAS in the structures considered give practically 
indistinguishable x-dependencies of the scattered wave amplitude in the grating. Therefore, the 
dependencies presented in figures 2(a)-(d) can equally be regarded as approximate and rigorous 
(within an accuracy of « 0.1%). 

Typical approximate dependencies of the normalized amplitude of the incident wave on the x- 
coordinate inside a narrow grating are presented in figures 3(a), (b) for the same structure as figures 
2(a), (b), i.e. with £2 = 5, e g = 5 x 10~ 3 , A = lum, 9iq = 45°, L = 10^m, and Ae3 = (i.e. there is no 
variation of the mean permittivity at the rear grating boundary). 

If the variation of the mean permittivity at the front boundary Aei i (i.e. ei < e 2 — figure 3(a)), the 
scattered wave does not carry the energy away from the grating, and energy conservation results in 
the same magnitudes of the amplitudes of the incident wave at the front (x = 0) and rear (x = L) 
boundaries (figure 3(a)). If Aei > (i-e. ei > e 2 — see figure 3(b)), then the scattered wave at x < is 
a propagating wave carrying energy away from the grating. As a result, the amplitude of the incident 
wave at the rear boundary is smaller than at the front boundary (figure 3(b)). 

If the mean permittivity varies at the rear boundary, i.e. Ae3 7^ and Aei = 0, then for small 




Figure 3: The x-dependencies of relative amplitudes |S2o(x)/Sio| of the incident wave inside the 
holographic grating with the same parameters as for figure 2 (£2 = 5, e g = 5 x 10~ 3 , L = lO^m, 
010 = 7r /4, A(vacuum) = ljum). (a, b) Aei = £\ — e 2 7^ and Ae3 = 63 — £2 = (variations of 
the mean permittivity only at the front boundary). Curves 1: Aei = 0, Ae3 = 0. Curves 2, 3, 4: 
(a) Aei = -10- 4 , -10~ 3 , -KT 1 , respectively; (b) Ae a = 10~ 4 , 10~ 3 , 10" 1 . (c, d) The approximate 
(thick solid curves) and rigorous (oscillating thin curves) dependencies of the relative incident wave 
amplitudes in the same grating but with (c) Ae3 = 0, Aei = 0.1 (curves 1) and Aei = _ 0-1 (curves 
2); (d) Aei = 0/ Ae 3 = 0.1 (curves 1) and Ae 3 = -0.1 (curves 2). 



variations of the mean permittivity (Aei ,3 <S £2) the approximate dependencies of the incident wave 
amplitude in narrow gratings are approximately the same as for the case with Ae3 = and Aei 7^ 
(figures 3(a), (b)). The smaller the grating width, the more accurate is this statement. Note, however, 
that this is not true for large variations of the mean permittivity Aei ,3 ~ e 2, for which reflection of 
the incident wave from the front boundary becomes noticeable. 

If Aei and/ or Ae3 are positive in a narrow grating with L < L c , then it is possible to choose optimal 
values of Aei ,3 <^ e 2 such that the amplitude of the incident wave at the rear boundary is next to zero. 
The smaller the grating width, the smaller the amplitude of the incident wave at the rear boundary 
can be. In this case, almost total conversion of the energy of the incident wave into energy of the 
scattered wave occurs. This is radically different from the conventional scattering in narrow gratings 
with small amplitude, where scattering is very inefficient and the scattered wave amplitude is small 
(as is the energy flow in it). 

If Ae3 = and Aei 7^ 0, then the rigorous x-dependencies of the amplitude of the zeroth diffracted 
order are practically the same as the corresponding approximate dependencies of the incident wave 
amplitude (figure 3(c)). The main distinctive feature of the rigorous dependencies is the weak fast 
oscillations with the period of the order of the wavelength (figure 3(c)). These oscillations are related 



to boundary scattering of the scattered wave at the grating interface x = L [10]. The wave resulting 
from this boundary scattering propagates in the negative x-direction as if it is a mirror reflected in- 
cident wave from this boundary. Therefore, mathematically the zeroth term in equation (2) includes 
both the incident wave and the wave caused by the boundary scattering [10]. Interference of these 
two waves results in a standing wave pattern represented by fast oscillations of amplitude of the 

1 /2 

zeroth diffracted order. The period of these oscillations, A/ (2e 2 cos 620), is in excellent agreement 
with the period of oscillations of the rigorous dependencies in figure 3(c). 

It is interesting that if variations of the mean permittivity occur at the rear grating boundary, i.e. 
Aeo, 7^ 0, then the amplitude of oscillations of the rigorous dependencies may be significantly larger 
(figure 3(d)). Increasing | Ae3 1 results in increasing amplitude of these oscillations. This is due to 
reflection of the incident wave from the rear grating boundary with Ae?, 7^ 0. Mathematically, the 
reflected wave is also included in the zeroth term in equation (2). Its amplitude increases with in- 
creasing I Ae3 1, which results in stronger oscillations of the resultant interference pattern. 
It follows that the difference between the approximate and rigorous curves in figure 3(d) can be 
reduced by including the reflected wave into the approximate analysis. Indeed, the total electric field 
resulting from the interference of the incident and reflected waves is given as 

[S20O) + s 2r exp(-2ik 20x x)} exp(i/c 2 ox* + ik Wy y - icvt), (9) 

where S2o(x) is the x-dependent amplitude of the incident wave in the grating, and S2,- is the am- 
plitude of the wave reflected from the rear grating boundary (as mentioned above, S2r is given by 
the Fresnel equations for reflection of the incident wave with the amplitude S2o(L) from the interface 
x = L). Wave (9) has a periodically varying amplitude (the expression in the square brackets). If we 
plot the x-dependencies of the magnitude of this amplitude for the gratings corresponding to curves 
1 and 2 in figure 3(d), the resultant approximate (oscillating) curves will be very close to the rigorous 
oscillating curves in figure 3(d). (Note that in the rigorous approach the term S2r exp( — 2i/c20xx) is 
formally included in the amplitude S^oix).) 

3.2 Wide gratings 

If the grating width L > L c , then the pattern of scattering changes significantly. Typical air- 
dependencies of the scattered wave amplitude inside the grating of critical width (L = L c f» 30um) 
are presented in figure 4. It can be seen that in this case the effect of varying mean permittivity on 
the scattered wave amplitude tends to be localized in the half of the grating that is adjacent to the 
boundary at which the variation occurs (figure 4). This tendency becomes much more obvious if we 
consider EAS in wide gratings, e.g. with L = 80um > L c « 30um — figure 5. 

If the variation of the mean permittivity occurs at the rear boundary (figures 4(c), (d) and 5(c), (d)), 
then the effect of this variation on the scattered wave amplitude is always localized within the dis- 
tance L c /2 near the rear boundary. Everywhere else in the grating, the scattered wave amplitude is 
hardly affected by the varying mean permittivity for small and large, positive and negative values of 
Ae3 (figures 4(c), (d) and 5(c), (d)). This is because the diffractional divergence can effectively spread 
the scattered wave only within distances ~ L c /2 [7-9,12], and thus the effects of any perturbations 
at the rear boundary (e.g. variation of the mean permittivity) can be felt only within these distances 
(figures 4(c), (d) and 5(c), (d)). 

If a variation of the mean permittivity occurs at the front boundary, i.e. Aei 7^ 0, then for negative 
values of Aei (figures 4(a), 5(a)) the situation is largely the same as in figures 4(c), (d) and 5(c), (d). 
However, if | Aei | is large (of the order of £2), then the scattered wave amplitude is affected (reduced) 
everywhere in the grating — see curves 5 in figures 4(a), 5(a). This is mainly because when Aei is large 
and negative, the amplitude of the incident TE wave transmitted through the boundary x = into 
the grating is noticeably less than the amplitude of the incident wave in the region x < 0. This must 
also result in a reduction of the scattered wave amplitude everywhere in the grating. 
If Aei > (figures 4(b), 5(b)), then increasing Aei results in a monotonous decrease of the scattered 
wave amplitude near the front boundary to about zero (similar to figures 4(a), 5(a)). However, in- 
side the grating (and at the rear boundary) the situation is different. The scattered wave amplitude 
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Figure 4: The x-dependencies of relative scattered wave amplitudes | S21 (z) / S\q | inside the holo- 
graphic grating with the same parameters as for figure 2, but with L m L c m 30,um. (a, b) Ae\ ^ 
and Ae 3 = 0, (c, d) Aei = and Ae 3 ^ 0. Curves 1: Aei = and Ae 3 = 0. Curves 2, 3, 4: (a) 
A £l = -10- 4 , -10- 3 , -10- 1 , respectively; (b) Ae a = lO" 4 , lO" 3 , 10 _1 ; (c) Ae 3 = -lO" 4 , -KT 3 , 
-10 -1 ; (d) Ae 3 = KT 4 , KT 3 , 10" 1 . Curves 5: (a) Aei = -4, i.e. e x = 1 (vacuum in front of the 
grating), and (b) Aei = 4.9 (in this case the angle of incidence 6\q = 45° is very close to the critical 
angle « 45.29° for the interface x = 0). 



decreases with increasing Aei (curves 1-3 in figures 4(b), 5(b)), goes through a minimum, increases 
back to approximately the same values as for gratings with e\ = ei = e 3 (see curves 1 and 4 in figures 
4(b), 5(b)), and then decreases again for large Aei (curves 5 in figures 4(b), 5(b)). 
The more complex behaviour of the curves in figures 4(a), (b) and 5(a), (b) compared to the curves 
in figures 4(c), (d) and 5(c), (d) is mainly related to the fact that if a perturbation occurs at the front 
boundary, then the effect of this perturbation can be spread into the grating not only by means of the 
diffractional divergence of the scattered wave, but also by means of the incident wave that propagates 
into the grating. At the same time, the effect of perturbations at the rear boundary can be spread into 
the grating only by means of the diffractional divergence of the scattered wave. Therefore the effect 
of varying mean permittivity at the front boundary may affect the wave amplitudes everywhere in 
the grating (figures 4(a), (b) and 5(a), (b)), whereas variations of the mean permittivity at the rear 
boundary may affect the wave amplitudes only within the distance « L c /2 from the rear boundary 
(figures 4(c), (d) and 5(c), (d)). 

The noticeable differences between figures 4(a), 5(a) and figures 4(b), 5(b) are due to the fact that if 
Aei > (figures 4(b), 5(b)), then the scattered wave in the region x < carries the energy away 
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Figure 5: Same as figure 4, but for grating width L = 80^im > L c « SOjUm. 

from the grating (see also the discussions of figure 2). The resultant energy losses from the grating 
can cause noticeable reduction of the scattered wave amplitude in the whole grating (see curves 1- 
3 in figures 4(b), 5(b)). However, if the variation of the mean permittivity Aei is sufficiently large, 
then the amplitude of the scattered wave at the front boundary is small (curves 4 in figures 4(b), 
5(b)). As a result, the energy losses drastically reduce, and the scattered wave amplitude inside the 
grating (not near the front boundary) increases — compare curves 2, 3 and 4 in figures 4(b), 5(b). For 
a particular value of e\ (in our examples this is e\ = e\ c = 10), the angle of incidence 9\q (if it is 
non-zero) becomes equal to the critical angle for the interface x = 0. In this case, the incident wave 
in the grating propagates parallel to this interface, i.e. 020 = 7t/2. On the other hand, if the angle of 
incidence in the grating 020 —> 7r /2, the amplitude of the scattered wave is strongly reduced [4]. This 
is the reason for the reduction of the scattered wave amplitudes represented by curves 5 in figures 
4(b) and 5(b). 

If the mean permittivity varies simultaneous at both the boundaries, then the effect of such variations 
is roughly the superposition of the effects of separate variations at each of the boundaries. This 
statement holds better for small variations of mean parameters. 

An important general feature of EAS can be seen from the approximate and rigorous analyses of 
scattering of bulk electromagnetic waves in narrow and wide gratings. Increasing or decreasing 
wavelength p times with the simultaneous increasing or decreasing grating width p times leaves the 
amplitudes of the scattered and incident waves unchanged (though scaled to different grating width). 
For example, EAS of bulk TE waves of A = 2/mn in gratings with L = 160 J um will be represented by 



exactly the same curves as in figures 5(a)-(d) with x/2 instead of x on the horizontal axis (scaling to 
the two times larger grating width). The other structural parameters are assumed to be the same as 
for figure 5. 

Similarly, p times increase or decrease of the dielectric permittivity in the structure (including grat- 
ing amplitude and step-like variations of the mean permittivity) with the simultaneous p 1 / 2 times 
decrease or increase of the grating width also leaves the wave amplitudes unchanged (though again 
scaled to the different grating width). 

4 EAS of guided modes 

As has been mentioned in the Introduction, the approximate method of analysis of EAS described is 
directly applicable to the case of EAS of guided optical modes of arbitrary polarization in gratings 
with constant mean structural parameters [2,5,6,9,12]. Here, we will see that this is also true for EAS 
of optical modes guided by a slab with varying mean thickness at the grating boundaries. In this 
case, the plane of figure 1 is the plane of the slab. One of the boundaries of the slab is periodically 
corrugated within the strip of width L (figure 1). Thus the thickness of the waveguide inside and 
outside the grating is given by the equation: 

( h, ( x<o, 

h(x ) = I h 2 + Z g f(x ) ,fort < x < L, (10) 
[ h, [ x > L, 

where h 2 is the mean thickness of the slab inside the grating region, hi and h$ are the slab thicknesses 
outside the grating, £, g is the grating amplitude (corrugation amplitude), /(xo) is a periodic function 
with the period A, max(|/(x )|) = 1/ and the average (over the period) value of f(x ) is equal to 
zero. Dissipation is neglected and all media in contact are isotropic. The variations of the mean slab 
thickness at the front and rear boundaries are given by Ah\ = hi — h 2 and A/23 = — h 2 (A/zi 3 can 
obviously be either positive or negative). 

As in the case of bulk electromagnetic waves in structures described by equation (1), the most inter- 
esting case of resonantly strong EAS of guided modes occurs at small grating amplitudes: 

|£gl/A«l- (11) 

In this case, the approximate theory (based on the two-wave approximation and the scalar theory of 
diffraction) is again expected to give highly accurate results (see [9,10,12]). 

Using speculations similar to those in section 2, one comes to the conclusion that equations (4) also 
describe EAS of guided modes in gratings with varying mean thickness at the grating boundaries. 
However, the coupling coefficients Tq and r\ are obviously different from those for bulk TE waves 
and are determined in the approximate theories of conventional Bragg scattering in corrugation grat- 
ings with grooves parallel to the grating boundaries [22-25]. The wave numbers for the incident 
and scattered modes (e.g. k\\, k 2 i, k 2 Q, etc.) are determined by the dispersion equation for the cor- 
responding slab modes. The effective dielectric permittivities corresponding to these modes can be 
introduced by the equations: £2ieff = k 2 ic/cu) 2 , £20eff = k 2 Qc/w) 2 , etc. 

Similarly to bulk electromagnetic waves experiencing reflection and refraction at the grating bound- 
aries with non-zero Aei 3, guided incident and scattered modes also interact with the step-like varia- 
tions of the mean thickness of the slab. For example, interaction of the incident mode with a step-like 
variation of h at the front boundary results in generation of other guided modes in the slab and bulk 
electromagnetic waves outside the slab [22-26]. In the approximate theory, however, all these waves 
do not contribute to scattering in the grating, except for one guided mode for which the Bragg con- 
dition is satisfied. Therefore, the situation is similar to the approximate theory of EAS of bulk elec- 
tromagnetic waves in section 2. For example, from consideration of the interaction of the incident 
mode with the front boundary (by means of the mode matching theory [22-26]), one can determine 
the amplitudes and the directions of propagation of the resulting bulk and guided waves, and then 



regard the guided mode with the amplitude S20O/ the angle of propagation 020/ and the wave vector 
k 2 o that satisfies the Bragg condition as the incident wave at x = in the grating (inside the grating, 
the amplitude of this wave is S2o(x)). 

It follows that equations (4)-(7) are directly applicable to the analysis of EAS of guided modes in 
corrugation gratings with step-like variations of the mean slab thickness. In this case, however, 
appropriate equations have to be used for the coupling coefficients, wave numbers, and it is assumed 
that S20, S21, Aq, and Bo are components of the electric or magnetic field in the incident and scattered 
modes. Note also that this approach is valid for both TE and TM guided modes. 
Interaction of the scattered wave with step-like variations of the mean thickness at the grating bound- 
aries also results in the generation of bulk and guided waves. This means that boundary conditions 
(8) (except for the first one) have to be modified in order to include these bulk and guided waves. 
This can again be done in the general case in the same way as in the mode matching theory [22-26]. 
However, manufacturing gratings with small amplitude usually results in small variations of the 
mean thickness of the slab. Adsorption and deposition of thin films on the slab surface (e.g. in the 
region x < 0) also results in only small variations of mean slab parameters (thickness or effective 
dielectric permittivity). Therefore, small variations of mean slab thickness are of most importance 
from the viewpoint of experimental observation and applications of EAS. 

The analysis of EAS in such structures does not require the general approaches of the mode matching 
theory and can be carried out by means of approximate boundary conditions similar to (8). This is 
because if Ah\ r 3 <: Ii2, then reflection and transformation of modes at such small stepwise variations 
of the mean slab thickness can be neglected. It can be seen that in this case, for TM modes, the 
boundary conditions at the grating boundaries x = and x = L can be written as 

Sio\ x =o = S2oo/E2izU=o = Eiiz|x=0/ (dE2iz/dx) x= o = (dEn z /dx) x=0/ 

E2iz|*=L = E3iz|*=L, (dE 2 iz/dx) x = L = (dE 3 i 2 /dx) x=L . (12) 

For TE modes, the z-components of the electric fields should be replaced by the z-components of the 
magnetic field. S20 is the amplitude of one of the components of the electric (or magnetic) field in 
the incident TM (or TE) mode inside the grating region (note that if Ahi <: h 2 , then S200 ~ S10 — the 
amplitude of the incident wave in the region x < 0, and 020 ~ $io)- The incident and scattered modes 
are assumed to be of arbitrary order. 

It can be seen that equations (12) are exactly the same as equations (8) with the only difference that 
the z-components of the electric or magnetic fields in equations (12) replace the total electric field in 
equations (8). 

When using the boundary conditions (12), it is convenient to use effective dielectric permittivities 
for the scattered modes: en e ff = k\\cj uo) 1 , £2ieff = kuc/tv) 2 , and £3i e ff = k^cjuS) 1 in front, inside 
and behind the grating, respectively. In this case, step-like variations of the mean slab thickness at 
the grating boundaries are related to the corresponding step-like variations of the mean effective 
permittivities. For example, for the front boundary: 

2k 1 \c 1 ( dfc 2 i\ A , 
Aen ef f = eneff - e 2 ieff ~ —^2— ) h h A/z i- ( 13 ) 

As mentioned above, this formula (and the similar one for the rear boundary) is correct only if 
A/*i,3 <C h 2 . Note that we do not have to consider step-like variations of the effective permittiv- 
ity for the incident wave. These variations can be neglected since for Ah\ ,3 <C I12 we have: S\o ~ S200 
and 020 ~ $10 (and similar for the rear boundary), i.e. interaction of the incident mode with the grat- 
ing boundaries can be neglected. (Note that, though being small, Aen e ff cannot be neglected owing 
to its significant effect on the diffractional divergence of the scattered wave.) 

For example, consider EAS of an incident TE zeroth (TE0) mode into the scattered TE0 mode in the 
structure: vacuum-GaAs slab (permittivity 12.25)-AlGaAs substrate (permittivity 10.24); the mean 
slab thickness within the grating region //2 = 6 x 10~ 5 cm, 020 ~ ^10 = tt/4 (for small values of Ahi <C 
\i2), the wavelength in vacuum A = 1.5 x 10~ 4 cm, the corrugation is assumed to be sinusoidal, 



i.e. f(xo) = sin(xo), and the effective permittivity £2ieff ~ 11.48. The period of the corrugation, 
determined by the Bragg condition, is A ps 0.579 J um. The wavelength of the guided TE mode in 
this structure is 0.42 J um, which is 92% of the wavelength (0.45jum) in the grating used for figures 2-5. 
Therefore, according to the general tendency discussed at the end of section 3, in order to get the same 
dependencies for the guided modes as in figures 2-5, we need to use the gratings of widths 0.92L 
(where L is the width in figures 2-5), and the step-like variations of the slab thickness determined by 
equation (13) with 

Aei 3 

Aei,3eff ~ -£21eff/ (14) 

where £2 = 5 and Aei,3 are the variations of the mean permittivity in figures 2-5. In this case, the 
corrugation amplitude can be chosen (£, g ps 2 x 10~ 6 cm) so that the x-dependencies of amplitudes 
of the scattered and incident slab modes in the grating are exactly the same as those given by figures 
2-5 with the replacement of x by x/0.92 on the horizontal axis (scaling to the new grating width). 

It can be seen from equations (13) and (14) that the variations of the mean dielectric permittivity in 
figures 2-5 Ae 1/3 = 1(T 5 , 1(T 4 , 1(T 3 , 1CT 1 correspond to Ah 1/3 w 1.3 x lCT 3 nm, 0.13 nm, 1.3 nm, 
130 nm in the considered guiding structure. Thus variations of the mean thickness of ~ 0.1 mono- 
layer may result in noticeable variations of the scattered wave amplitude in the structure — see curve 
2 in figures 2(b), (d). Obviously, on the one hand, this is a significant complication to the experimental 
observation and application of EAS of guided modes. On the other hand, this demonstrates a very 
high sensitivity of EAS to conditions on slab interfaces, which may open up excellent opportunities 
for the design of highly sensitive optical sensors (such as adsorption sensors, sensors of dielectric 
permittivity, etc.). 

If the high sensitivity of EAS to varying mean structural parameters is not desirable, then one should 
use wide gratings, where the effect of varying mean parameters at the grating boundaries is notice- 
able only within a half of the critical distance — see figures 4 and 5. Also, use of larger grating am- 
plitudes results in reduction of the sensitivity of EAS to varying mean parameters. Another option 
is to use grazing-angle scattering (GAS) [12] rather than EAS, where the amplitude of the scattered 
wave may be especially large at the middle of the grating but not at its boundaries. There is also a 
possibility of compensating (at least partial) for varying mean parameters by choosing the angle of 
scattering so that the scattered wave outside the grating propagates parallel to the grating [27]. 

Another option for the reduction of the high sensitivity of EAS for guided modes is to use slabs with 
larger thickness h, and permittivity that is closer to the permittivities of the surrounding media. For 
example, in the structure with a polymer slab (dielectric permittivity 2.56) on a silica substrate (dielec- 
tric permittivity 2.22), h = Ifim, A = l^jtmi, L = 20jum, and grating amplitude £ ? = 3.3 x 10~ 5 cm we 
have: £2i e ff = 2.49, and the equivalent structure for bulk waves (having the same dependencies of the 
wave amplitudes) must have £2 = £2i e ff/ same A and L, and e g = 7 x 10~ 3 . In this case, A£i = 10~ 4 in 
the equivalent structure corresponds to a variation in the mean thickness A/zi ps 1.8 nm. These varia- 
tions of the permittivity or thickness result in 15% reduction of the scattered wave amplitude from 
~ 8S200 to ~ 6.8S200 a t the front grating boundary x = 0. Recall that a similar 15% variation of the 
scattered wave amplitude in the GaAs waveguide with L = 9.2jUm, h = 0.6,um, and ^ = 2x 10~ 6 cm 
is achieved at a much smaller value of Ah\ pa 1.3 x 10~ 2 nm — compare curves 1 and 2 in figure 2(b). 
Finally, it is important to note that in the case of EAS of surface waves (e.g. surface acoustic waves, or 
surface electromagnetic waves) in corrugation gratings, the problem with varying mean parameters 
hardly exists. This is because in this case the step-like variation of the mean level of the substrate 
surface inside and outside the grating does not affect the length of the wave vector (or the wave- 
length) of the surface waves. Only if, for example, outside the grating we have additional layers on 
the surface, may the mean propagation parameters change and the above-mentioned effects occur. 
Thus a surface wave sensor based on EAS in periodic groove arrays is especially easy to design. 



5 Conclusions 



This paper has demonstrated yet another unique feature of scattering of optical waves in extremely 
asymmetrical geometry — the unusual sensitivity of the incident and scattered wave amplitudes to 
small stepwise variations of mean structural parameters (e.g. mean dielectric permittivity or mean 
waveguide thickness). This effect has no analogies in conventional Bragg scattering. It has been 
explained by the high sensitivity of diffractional divergence to small variations of mean structural 
parameters across an optical beam. 

Two distinct typical patterns of scattering have been described for gratings that are narrower and 
wider than a grating of critical width L c [7-9] (half of L c is equal to the distance within which the 
scattered wave can be spread across the grating by means of the diffractional divergence, before it 
is re-scattered by the grating [7-9]). In particular, it has been shown that in narrow gratings (with 
L < L c ), varying mean permittivity at either of the boundaries strongly affects the scattered wave 
amplitude everywhere in the structure, whereas in wide gratings (with L > L c ), the effect of vary- 
ing mean parameters is usually (but not always) significant only within the distance L c /2 from the 
boundaries. 

The analysis has been carried out by means of approximate and rigorous methods, showing a very 
good agreement between them for the most interesting cases of strong EAS in gratings with small 
amplitude. 

The analysis of EAS of guided modes has revealed especially high sensitivity of scattered wave am- 
plitudes to small variations of the mean slab thickness. This is mainly due to strong dispersion of 
guided modes, i.e. strong dependence of their wave numbers on slab thickness. This sensitivity may, 
on the one hand, present a complication for experimental observation of EAS, and on the other hand, 
be very useful for the application of EAS in the design of highly sensitive sensors and measurement 
techniques. Several options for dealing with this possible experimental complication are discussed. 
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